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Abstract 

The classical Bloch equations of spin magnetization transport is extended to fully time-dependent 

d ' and highly-nonlinear nonequilibrium spin magnetization quantum distribution function transport 

(SMQDFT) equations. The relevant variables are the spinor correlation functions which separate 

. into charge and spin magnetization distributions that becomes highly coupled in SMQDFT equa- 

>. 

tions. The leading terms consist of the Boltzmann kinetic equation with spin-orbit coupling in a 

(N 

CN ■ magnetic field together with spin-dependent scattering terms which contribute to the torque. These 

■ do not have analogue within the classical relaxation-dephasing picture, but are inherently quantum 

many-body effects. These should incorporate the spatio-temporal-dependent phase-space dynam- 
ics of Elliot- Yafet and D'yakonov-Perel scatterings. The resulting SMQDFT equations should serve 
^ ■ as a theoretical foundation for computational spintronic and nanomagnetic device applications, in 

5^ . ultrafast-switching-speed/low-power performance and reliability analyses. 
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The ultrafast-switching-speed and power-dissipation performance analyses of nano- 
electronic devices has revealed the inadequacy of the classical Boltzmann transport equa- 
tion (BTE) and ushered its extension of the to fully time-dependent and highly-nonlinear 
nonequilibrium quantum distribution function (QDF) transport equations for charge carri- 
ers. This has been achieved through the use of non-equilibrium Green's function, obtained 
either by the time-contour quantum field formulation of Schwinger[l|, KeldyshQ], and 
Kadanoff and Baym 3| , or by the real-time quantum superfield formulation of Buot cou- 
pled with his lattice Weyl (LW) transformation technique {4- 3]- The extension of BTE to 
QDF transport equation has proved to be highly crucial in discovering autonomous current 
oscillations in resonant tunneling devices through numerical simulations, and in resolving 
controversial issues concerning highly- nonlinear and bistable current- voltage characteristics 
found in the experiments jol. [lo|. Indeed, the QDF approach in phase space has so far been 
the most successful technique in the time- dependent analyses of open and active nanosystem 
and nanodevices[lo|. 

To the authors' knowledge, the extension of the well-known classical Bloch equation 
for spin transport, which is the analogue of the classical BTE for particle transport, to 
fully time-dependent and highly-nonlinear nonequilibrium QDF transport equations for the 
magnetization has not been reported. With the exploding surge of interest on spintronics 



and nanomagnetics 



111 ], there is an urgent need for this fully quantum transport extension of 



the classical Bloch equation to guide the numerical simulation of the speed-power switching 
performance of realistic spin nanostructures/transistors. The well-known classical Bloch 
equations for spin systems is usually given in the diffusive regime 12] as 
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where 7 is the gyromagnetic ratio, Ti is the spin-relaxation characteristic time, and T2 is 
the spin-dephasing characteristic time. 

The purpose of this paper is to extend this classical Bloch equations to a fully time- 
dependent and highly-nonlinear nonequilibrium SMQDFT equations on the same level as in 
the extension of the classical BTE to full QDF transport equation based on nonequilibrium 



Green's function technique [l-l6|. In this paper we will first treat spin transport within a 
single energy band. The multiband spin transport, where pseudo-spin plays a role, will be 
reported in a separate arxiv.org online paper. 

Our starting point is the general quantum transport expressions for fermions, as ob- 
tained from the real-time quantum superfield theoretical technique js], Q]. Upon dropping 
the Cooper pairing terms in those expressions which do not concern us in this paper (their 
corresponding transport equations j^, O] are important in nonequilibrium superconductiv- 
ity), the nonequilibrium Green's function transport equations reduce to the following known 
expressions P 
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+ [S^G^ - G^S"] + [S^G'^ - G"S^] . (2) 

The quantum distribution-function transport equations in (p, q) = (p, q, E, t)-space are 
obtained by applying the LW transformation |3] (although continuum approximation is em- 
ployed in this paper, this is not essential and we use the word "lattice" when referring to 

n n 

solid-state problems) p, |6[. The LW transform a{p,q) of any operator A is defined by the 
following relation 

a{p,q) = y"rfz;exp (^^^ = W A(v,r,q,t), (3) 

where the matrix element in the integrand is evaluated between two Wannier functions, and 
W indicates the operation of taking the LW transform. This transformation enables the 
numerical simulation of real time-dependent dynamical open systems by transforming from 

two-point correlation functions to numerically manageable and physically meaningful local 

I 

functions in phase space [8f|. The well-known electron Wigner distribution function (WDF), 
fw (p, q, t), is obtained from the correlation function —ihG"^ by performing integration over 
the energy variable 

U (P, q,t) = ^JdE {-ihG< (p, E, q, t)) . (4) 

In the absence of particle pairing of superconductivity to obtain Eq. ([2]), we may also 
write the general transport equation,written for '< '-quantities as [here {A^B} and [A^B] 



means anticommutator and commutator, respectively, of A and B] 



+ [S<,ReG'1-^{r,G'<} + ^{S<,A}, (5) 



where the single-particle Hamiltonian, 

= Ho + S" + Re S 

with E*^ being the 'singular part' (delta function in time) of the self-energy. 

'Ho = Eip) + Viq), 

where E (p) is the electron dispersion relation or energy-band function and V [q) is the 
self-consistent single-particle potential function. In obtaining Eq. from Eq. we 
made use of the relations, (G"^)^ = C, and S^'''(l,2) = — S^(2,l), and the follow- 
ing relations, iA{l,2) = -2zlmG^ = - (G> (1, 2) - G< (1, 2)), ir(l,2) = -2iImS^ = 
-(S>(1,2)-S<(1,2)). 

The second term in Eq. ([S]) is a commutator and also describes evolution in phase- 
space although only serving to complicate, by virtue of the scattering term, S^, the kinetics 
of particle motion described by the first term, a sort of complex interference phenomena 
similar to zitterbewegung in Dirac quantum mechanics 13|] . Although often neglected in 



considering quantum transport of particles as this is expected to make a small corrections 
to the particle current, we will include the [S^,ReC] term in our present considerations 
for completeness. Note that in the gradient expansion of Eq. (E]), with the [E<, Re C] term 
neglected, the leading terms immediately give the classical Boltzmann transport equation 



for spinless systems IJ] . However, even this is no longer quite true for spin system as we 
shall see later, since the spinor scattering terms contribute to the torque in the system. 

The third term of Eq. (|5]), which is proportional to G^, is the scattering- out term, and the 
last term which is proportional to the spectral density is the scattering-in term, in analogy 
with the collision terms in the classical Boltzmann equation. The collision terms can also 
be rewritten using the identity, 

-\ {r, G<} + ^ {S<, A} = -\ {S<, G>} + \ {S>, G<} , 



which also leads to iS^ and — iS*^ terms, upon performing the W operation, as scattering 
out and scatterine:-in 



The equation for G"" and = G"'^ are simpler and are given as, G ^G^ = 0, G"" [G ^ 
0, where 



G-'G' (1,2) = ih ( ^ + ^ ] G' (1,2) - U + G' 

Oti Ot2 
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The above equations are then cast in the language of phase-space QDF transport equa- 
tions, similar to that of the classical Boltzmann equation by the use of Buot lattice Weyl- 
Wigner formulation of the band dynamics of electrons in a solid In the continuum 
approximation, the first thing to do is to recast the two-point space and time arguments in 
all the transport equations as follows: ti = t— t2 = t+'^, (\i = ^2 = We also 

need to define a (3 + l)-dimensional canonical variables, namely, p = (p, —E) , g = (q, t) , 
which allows us to write v = (v, r). The transformation of Eq. (j3]) is then applied. 

For simplicity in spin quantum transport equations, we consider only a single energy band 
having spin indices, I and Incorporating the J, and f variables in Eqs. (E]) and ([5]), we 
obtain exactly four transport equations for G^^ (1, 2), (1)2), G^^ (1, 2), and G^^ (1,2). In 
spintronics, we are interested in the time-dependent evolution of the magnetic polarization, 
Sz = (G^^ — G|^) , as this is transported across the device. In order to achieve this, we 
need to perform linearly-independent combinations of the components of the spinor nonequi- 
librium Green's function consisting of the set: {G^^ (1, 2) , G^^ (1,2), G^^ (1, 2) , G^ (1, 2)}. 
The resulting new set will become the relevant and independent correlation functions more 
pertinent to the spin-transport problem, which consists of the separation into charge distri- 
bution and spin magnetization vector. 

The separation in terms of a scalar and vector representing the total charge and spin- 
vector correlation functions, naturally occurs as expansion coefficient in terms of the Pauli 
matrices. 
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where / is the 2x2 identity matrix, and 




> 



(7) 



where we drop the '<' superscript in the spin correlation functions, Sj. In the frame of 
reference where the z direction is fixed by the magnetic field, the and Sy evolution 
equation describe dephasing mechanisms in the x-y plane of the 'Bloch sphere'. The total 
charge of the system, p {p, q, t) is represented by the scalar So i.e., from Eq. (|4]) we have. 



where So {p, E, q, t) is the LW of correlation So- Thus Eq. performs the separation of the 
charge and spin correlation functions as the relevant variables of any fermion system with 
spin degree of freedom. 

For slowly varying system where the total charge vary slowly with space and time, per- 
haps in some spin transfer torque systems, we can reduce the number of important spin 
correlation functions in the spin quantum transport problem by making an approximation 
to S'o[15]. Assuming the conservation of particles possessing spin degree of freedom within 
semiconductor channel of, say, a spintronic transfer torque transistor, we must have the fol- 
lowing relation for the 'tf and 'J,^' rates of change, ^^^^gf'^'*'' = — ^^^ii^^^ where p^f (p, q, t) 
and pii (p, q, t) are the LW transforms of ihG^^ and ihG^^, respectively. 

Generally however. So will be varying in both space and time, i.e., a two-point nonequi- 
librium correlation function in space and time, describing the varying charge within the 
spintronic device due to depletion and accumulation in various regions within the device, 
such as in spin transfer torque devices employing thin barriers between the channel and 
ferromagnetic leads. 

We can also write the spinor self-energies used in Eq. ([2]) asQ, Q into spin-independent 
scalar [denoted by overscript bar symbol] and vector quantities. 
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which we write exphcitly as, 
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Likewise, the spinor quantities used in Eq. ([5]), namely, the spin-dependent scattering-out 
rterm is cast in the form, 



= (imE'^J + ImS^ -a) , (9) 

so that 

r = -ImS^ 

7 = -ImS''. (10) 
Similarly, we have for the spin- dependent spectral function. A, in the scattering-in term as 

Aap 1 



-(^Ai + A-a')=-lmS:^ 



ilmS^'I + lmS'' - a) , (11) 



2 

so that 

A = -Im^", 

^=-ImS^^ (12) 

In the effective-mass approximation, we may take the corresponding single-particle Hamil- 
tonian Tig for spin systems in a magnetic field with spin-orbit coupling as 

•H.= ( ^ - eV {q) + ^Ref:Ai + H^o (v, q) + ^ (Rc S^' ■ a - gefffieB ■ a) , (13) 
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where V = P + (g), electric charge, e = \e\, V (q) is the self-consistent applied potential, 
and B = Vg x A{q). The second term, Hsnjv, , is the spin-orbit coupling which may 
consists of Dresselhaus and/or Rashba term 18| as well as transistor-gate and strain-induced 



spin-orbit coupling. The term Re S*^' ■ a comes from the spin vector of the spin-dependent 
part of Re S^^ in the single particle Hamiltonian. Its presence basically modifies the effective 
magnetic field due to external fields and spin-orbit coupling. 

We have in general. 



qx V 
2hm*c^ 



E 



which may acquire definite form in lower dimensional system based on symmetry considerations 
In what follows, we write the most general form of the single-particle Hamiltonian, T-Lsi ^ 

where h {P,Q^ = (^gs^tAf ' ^^^S^) . Therefore 

Hs = \{h i-a-B), (14) 

where 

and for electron-spin angular momentum, the g-factor, Qs — 2. Here the effective B^ff (^P, Q 
includes the effects of transistor-gate or strain induced spin-orbit coupling, and applied mag- 
netic fields. H (^P, consists of the energy-band structure and spin-independent confining 
potential due to the device heterostructure coupled with the applied bias at the terminals. 
The spin-orbit coupling in B^fj ^-P, also depends very much on the electric field of the 
self-consistent potential inside the device. This in turn depends on the total electric charge, 
represented by So, through the Poisson equation. Although we could include Re S*" in the 
definition of B above, we treat this separately to emphasize the possible role of scattering 
terms in contributing to the torque in the system. 

The resulting equations obeyed by the spin correlation functions defined by Eq. ([7]) in 
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accordance with Eq. ([2]) can be written in vector notation as, 
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and for the scalar correlation function representing the total charge , 

= \ [H, s<] + \ [f.^s< - s<f.^] + \ [t.<s: - s:t.<] 



+ 



B-S< -S< -B 



S^-S 



r --:< 



(16) 



(17) 



For self-consistency, one must also be solved for potential, or electric field, E, using the 
Poisson equation. 



P {q, t) - po 



where $ = V in Eq. (fT3|) . po is the positive background charge and p{q,t) is given by Eq. 



We also have the equivalent equation for 5"* 



_ 1 

~ 2 

i 

+ 2 
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B + ReE'] X S< - S< X (B + Re 

1 



j:<,ReS' 



- ^ jlmS" X 5< + S< X ImS'l + ^ |S< x Ini^" + Im S^'^ x S<} 
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H< X Re^" - Re5" x H< 
1 
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ReS'",5< 



S<,Re5: 



--^S<,Im5;^, (18) 



which can be reduced to that which one can derive from Eq. ([5]), by using the relations in 
Eqs. ([TO]) and (in]). The resuh is 
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dti dt2 



_ 1 

~ 2 

i 

+ 2 



+ 2i7 y S'-+ S'=X7 



< c'< 



i3 + ReSM X 5< - S< X m + ReS 



1 



+ 



+ 



X Re 5'^ - Re S'' x 



i3+ReSM ,5, 



+ 



Re 5! 
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-\{l.s<}-'-{'^<^ms: 



(19) 



Note that our definition of i3, Eq. ( ITSl) . has the unit of energy, just hke the self-energies. 
The equivalent expression derived from Eq. ^ in terms of the magnetization vector, 
= 9^2^"^ =9sfJ'B (the ^(-factor for spin, gs ~ 2) is. 



S<,ReM' 



2 

+ 2 



i3 + ReSM X 7W - 7W X n3 + ReS" 



+ 



+ 



S< X Re>l" - Re>l" x S< 

1 



+ 



H<,ReM: 
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(20) 



where we defined, 



First, let us evaluate Eqs. (120|) by assuming constant magnetic field B and in the absence 
of scatterings. We obtain 



d_ 

di 



M 
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(21) 
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which is the balhstic version of the form of Eq. ([T]) . 

In carrying out the LW transform, the following relations between vector cross products 
and commutation or anti-commutation relations are useful, 

Ax B - B X A = i.eijk {A^^Bk] , 
Ax B + B X A = iiiijk [Aj, Bk] , 

where /j is the unit dyadic symmetric tensor or idemfactor, and iijk is the anti-symmetric 
unit tensor. For the purpose of making gradient expansion to make contact with classical 
expressions, we give here the following for convenience the LW transform of a commutator 
[A, B] and an anticommutator {A, B} in terms of Poisson bracket differential operator. A, 
as 

[A B] {p, q) = cos A [a {p, q) b {p, q)-h {p, q) a {p, q)] 

- isinA {a {p,q)b{p,q) + b {p, q) a {p, q)} , (22) 
{A, B} (p, q) = cos A {a {p, q) b {p, q) + b {p, q) a {p, q)} 

- z sin A [a {p, q) b {p, q) - b {p, q) a {p, q)] , (23) 

where A = | ■ ^ — ■ • Thus if the LW transforms are not matrices, then to 
lowest order, we have 

[AB\[p,q}--i-\^-^ - - ^J, (24) 

{A, B} {p, q) = {a (p, q) b (p, q) + b (p, q) a (p, q)} = 2a (p, q) b {p, q) . (25) 

If we compare Eqs. (IT6l)- (fT9l) with the particle quantum transport equation of Eq. ([5]), 
without the zitterbewegung [S^, Re C] term, we see that the leading terms of the first line of 
the RHS of Eq. f|T9|) . corresponds to the classical BTE. The next three lines give the torque, 
due to spin-orbit coupling (included in the effective magnetic field), and terms coming from 
the spin-dependent relaxation and dephasing scattering mechanisms. The last line explicitly 
involves coupling to the total charge. 

Carrying out the W operation in Eq. (|20|) . we find that the leading gradient terms in the 
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RHS reads 
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dp dq dq dp 

where all quantities above are LW transforms, i.e., they are functions defined in phase space 
(p, t) = {p, —E, q, t). In the first line of Eq. fl26|) and the first term of the second line arise 
from the relation, 



'd{H + ReJ:'^) dM d{H + ReE') dM 



dp 



dq 



dq 



m 



p 



dp 
V X B 



dReJ:'-dM dReE"- dM 



dp dq dq dp ^ ^ ^ 

where the RHS arise from the expression given by the first terms of Eg . ([1311 . We have used 
the gauge: A{q) = x B is the uniform external magnetic field [19]. 

The first line in the RHS of Eq. ( l26l) is the Boltzmann equation in a magnetic field with 
effects of spin-orbit coupling. The second line is the contribution to the Boltzman equa- 
tion if the effects of Re and ReA^^ are included. The rest of this equation incorporates 
spin-dependent scattering, ReS'', which directly contribute to the torque. The precession 
of the magnetic moment, Ai, is also correlated with the torque exerted by 7, and similarly 
for ^ scattering-in^ precession (fluctuation) of A due to torque exerted by E. The last two 
lines describe the coupling to the total charge. Thus, spin-dependent scattering terms are 
capable of describing the spatio-temporal-dependent scattering dynamics of various mecha- 
nisms of spin relaxations and precession, namely, Elliott- Yafet (spin-dependent impurities 
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and phonon scatterings), D'yakonov-Perel' (due to electron-electron scatterings, with elec- 
trons executing random- walk in momentum space, resulting in the fluctuations of magnitude 
and direction of spin precession axis). The Bir-Aronov-Pikus scattering (due to the interac- 
tion with holes in valence band, serving as a spin sink), and the hyperfine interactions are 
not relevant in our present treatment. 

Until now the treatments of spin transport are mainly focused on scattering mechanisms, 
based on quasi-classical steady-state conditions and/or where nonuniformity in real space is 
often not considered, common in micromagnetics. On the other hand for spintronic transistor 
devices, strong inhomogeneity in real space is of prime considerations. The calculations of 
the self energy S often have the goal of obtaining the relaxation times to be used in classical 
kinetic equations. 

A calculation where the self-energy spinor due to many-body effects of spin-orbit coupling 
is treated has been given by Rajagopal|21[. A first principle ab initio calculation of the self- 



energies has been given by Zhukov and Chulkov|22l] for the Elliot- Yafet mechanism in Al, Cu, 
Au, Nb, and Ta, and by Mower, Vignale, and Tokatly 23| for D'yakonov-Perel spin relaxation 
mechanism in photo-excited electron-hole liquid in intrinsic semiconductors exhibiting spin- 
split band structure. 

Our derivation of the nonequilibrium spin magnetization QDF transport equations is 
based on the separation of charge and spin as the relevant variables. The resulting SMQDFT 
equations should provide the fundamental basis for carrying out the numerical simulation of 
the switching or time-dependent performance analyses of spintronic devices, such as those 
that make use of thin insulating layers between conducting metal layers [24J]. Potential bar- 
riers of up to 1 eV exist in promising spintronic devices, exemplified by the spin valve and 



magnetic tunnel transistors 



25 
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